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Abstract: We analyze the antimatter yield of supersymmetric (SUSY) models with large 



f^ , neutralino annihilation cross sections. We introduce three benchmark scenarios, respec- 



tively featuring bino, wino and higgsino-like lightest neutralinos, and we study in detail 
the resulting antimatter spectral features. We carry out a systematic and transparent 



p i' comparison between current and future prospects for direct detection, neutrino telescopes 

O ; and antimatter searches. We demonstrate that often, in the models we consider, antimat- 

^ ■ ter searches are the only detection channel which already constrains the SUSY parameter 

space. Particularly large antiprotons fluxes are expected for wino-like lightest neutralinos, 
while significant antideuteron fluxes result from resonantly annihilating binos. We intro- 
rS \ duce a simple and general recipe which allows to assess the visibility of a given SUSY model 

C^ I at future antimatter search facilities. We provide evidence that upcoming space-based ex- 

periments, like PAMELA or AMS, are going to be, in many cases, the unique open road 
towards dark matter discovery. 
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1. Introduction 

Weakly interacting massive particles (WIMPs) are among the leading candidates for cold 
dark matter in the Universe: on one side, the idea is appealing since dark matter appears 
just as another of the thermal leftovers from the early Universe; at the same time, from the 
particle physics point of view, numerous models naturally embed such kind of particles. 
The search for WIMPs has therefore been one of the hot topics in experimental physics over 
the last decade or so, and the field is now acquiring new perspectives with next-generation 
experiments being at the research and development stage, or, in some cases, even closer to 
the data taking phase. The target is either to detect WIMPs directly, measuring the recoil 
energy from WIMPs elastic scattering on nuclei, or to search for the yields of WIMP pair 
annihilations, such as neutrinos from the Earth/Sun or cosmic rays from annihilations in the 
galactic halo. Complementarity between direct detection and indirect detection methods 
has been repeatedly stressed (for reviews on this topic, see, e.g., [|l|, |2|). Indeed, the idea of 
WIMP dark matter refers to a broad framework rather than to any specific model; in such 
a framework, an overall (weak) coupling of dark matter to ordinary matter is guaranteed. 



but details on how such couphng is reahzed may be sensibly different from one model to 
another: One of the most pressing issues from an experimentalist perspective, the question 
regarding how competitive a given detection technique is compared to the others, can be 
addressed just on model dependent grounds. 

The most popular WIMP dark matter candidate is the lightest neutralino in supersym- 
metric extensions of the standard model of particle physics. The idea of supersymmetry 
has strong theoretical motivations, but it has, unfortunately, no experimental confirmation 
so far. Recipes on how to introduce a super symmetric model (and hence the correspond- 
ing dark matter candidate) have rather weak constraints; to some extent, this freedom is 
even too large, so that, for practical purposes, it is customary to focus on minimal models, 
i.e. models with minimal particle content and a small number of free parameters. The 
benchmark scenario most often implemented is the mSUGRA setup, i.e. the supergravity 
inspired model with the smallest possible number of free parameters, namely only four 
parameters and one sign, partly defined at the grand unification (GUT) scale and partly 
at the low energy (electroweak) scale. Referring, within this framework, to models with a 
thermal relic density (in a standard setup) matching the cold dark matter term as deduced 
from latest high precision cosmological measurements Q, prospects for direct detection 
look more promising than those for indirect detection methods (see, e.g., the recent anal- 
ysis in ref. [^). To some extent, a similar conclusion is reached also for thermal relic 
models in the context of minimal supersymmetric standard model (MSSM) setups defined 
by few universal parameter chosen at a low energy scale (such as the seven parameter 
model implemented, e.g., in ref. [^). There is often, in the literature, the tendency to 
extrapolate on these results, and assume that, as a generic trend for supersymmetric dark 
matter candidates, direct detection should have top priority. 

We give here a few counter examples, still in the context of minimal supersymmetric 
setups, which allow for a clearer presentation of results and a more transparent comparison 
between different signals (for related studies within various supersymmetry breaking sce- 
narios see, e.g., B). We will enlighten, in particular, the role of indirect detection through 
the search for exotic antimatter components in cosmic rays. This technique was proposed 
soon after the formulation of the idea of WIMP dark matter, the links to the framework 
definition being evident: Pair annihilation is the mechanism which sets the thermal relic 
abundance of WIMPs; although the density of WIMPs in dark matter halos today is much 
smaller than in the early Universe environment, there is still a finite probability for WIMPs 
in the Galactic halo to annihilate in pairs. In these annihilations the same amount of matter 
and antimatter is produced; while the matter component is likely to be very subdominant 
compared to standard astrophysical sources, there seems to be no standard primary source 
of antimatter, with the bulk of the (scarce) antimatter component in cosmic rays which is 
likely to be of secondary origin, i.e. generated in the interaction of primary cosmic rays 
(mainly protons) with the interstellar medium (mainly hydrogen and helium). The goal 
is then to identify, through their peculiar spectral features, the WIMP-induced antimat- 
ter fluxes, or at least to exclude those dark matter candidates which would overproduce 
antimatter compared to the relatively low background term. 

Implementing models for the dark matter halo with self-consistent density and velocity 



distributions, we will produce fully consistent comparisons of indirect detection versus 
direct detection, examining both current exclusion limits and future detection prospects. 
Our analysis will be focussed on a few sample particle physics scenarios, but to facilitate 
the extension to alternative frameworks, we will introduce a simple and general recipe to 
address detectability of a WIMP model in future antimatter cosmic ray measurements. 

The outline of the paper is as follows. After introducing the particle physics setup in 
sec. |2|, we provide, in sec. ^, details on the computation of antimatter fluxes from neutralino 
annihilations, and describe the relevant spectral features. Sec. |^is devoted to a comparison 
among current exclusion limits from different search strategies, while in sec. y we present 
a novel technique for the evaluation of future discrimination capabilities at antimatter 
searches, and confront the projected experimental sensitivities at direct and indirect future 
dark matter detection experiments. The question of the dependence of antimatter yields 
on the halo model is treated in sec. |^ Finally, sec. ^ concludes. 

2. The particle physics framework 

We suppose that our cold dark matter candidate is the lightest neutralino x^ as defined in 
a MSSM setup, i.e. the lightest mass eigenstate from the superposition of the two neutral 
gaugino and the two neutral Higgsino fields: 

X? = NuB + NuW^ + Ni3H^ + NuH^ . (2.1) 

The bino, wino and higgsino fractions, which we define, respectively, as |iVii|, I-/V12I, 
andy^Nf^^ + N^^, are obtained by diagonalizing the neutralino mass matrix, and depend 
mainly on the relative hierarchy between the bino, the wino and the higgsino mass parame- 
ters, i.e., respectively, Mi, M2 and fi. The WIMP relic abundance and detection prospects 
are mainly set by the lightest neutralino mass and composition, but depend in part also on 
the structure of the Higgs and sfermion sectors of the theory. We will select a few specific 
configurations, and to clarify how this choice is made, we need to briefiy review which role 
is played by the relic density and by other constraints on the model. 

2.1 Relic abundance, accelerator constraints and WIMP models 

Let's start our discussion in the context of the popular minimal Supergravity (mSUGRA) 
framework [^. In this scenario, the lightest neutralino is 6ino-like over most of the param- 
eter space, as a consequence of gaugino soft SUSY breaking mass unification at the GUT 
scale; the only exception is a rather small parameter space region, often dubbed focus-point 
or Hyperbolic Branch region Q, close to the region where electroweak symmetry breaking 
(EWSB) conditions can no longer be fulfilled, where a low value of the // parameter is 
achieved through large values of the scalar mass parameter tuq. 

Binos have rather low pair annihilation rates (mainly into fermion-antifermion pairs, 
much smaller than higgsinos or winos annihilation rates into weak bosons), getting lower 
and lower for increasing bino masses. As the WIMP relic abundance scales with the inverse 
of the annihilation rate (the higher the rate, the longer the WIMP is kept in thermal equi- 
librium, the smaller the WIMP density at thermal decoupling and hence today), it follows 



that in most of the mSUGRA parameter space, the Ughtest neutrahno rehc abundance is 
larger than the upper WMAP bound on r^cDM^^- One should then restrict to light binos, 
in the so-called bulk region at small values of tjiq and of the soft super symmetry breaking 
fermionic mass parameter mi/2 (the one through which Mi, as well as M2 — 2Mi, is fixed 
at the low energy scale); however, the bulk region is almost entirely excluded by accelerator 
constrains since, at low values of tan (5, the ratio of the vacuum expectation values of the 
two neutral components of the SU(2) Higgs doublets, the LEP constraint on the lightest 
CP-even Higgs boson mass nifi is violated, while, at large tan/3, stringent bounds stem 
from the inclusive b ^ s'j branching ratio (for reference, throughout this analysis we will 
assume as a bound on m/j the standard model Higgs bound of 114.1 GeV y], and an allowed 
range 2.0 x 10~^ < BR[i? —>■ Xg 7] < 4.6 x 10~^, obtained adding a theoretical uncertainty 
of ±0.5 X 10~^ to the experimental value quoted by the Particle Data Group 2002 [^], with 
an estimate of this process including the complete next-to-leading order (NLO) correction 
for the standard model contribution [0, |l2| and the dominant NLO corrections for the 
SUSY term ||, Q). 

The loopholes for bino dark matter in this scenario are found in two regimes: 1) the tiny 
coannihilation branches, i.e. stripes of the parameter space where the lightest neutralino 
is just slightly lighter than the next-to- lightest SUSY particle, i.e., in this scenario, a stau 



at low mo [^, |l^ or a stop for large trilinear couplings |17|, which drives the thermal 
decoupling in the early Universe; 2) the funnel region, where annihilations into fermions 
through s-channel exchanges of the CP-odd A (or of the heaviest CP-even H^) neutral 
Higgs boson are enhanced by the resonance condition m^ ~ inAf^ (or m^ ~ m^Q/2). In 
both regimes, and in the lower neutralino mass end, direct detection rates are well below 
current sensitivities, but above the projected sensitivity of future large-mass experiments. 
Unfortunately this does not follow from crossing symmetry, i.e. from some diagrammatic 
structure giving large contributions to both the annihilation rate and the scattering rate 
on light quarks. The scalar superpartners of light quarks are very heavy, hence scattering 
rates are dominated by t-channel exchanges of two CP-even Higgs; these contributions are 
guaranteed to be large in this framework, as the higgsino fraction is never forced to zero 
(the coupling of neutralinos to a Higgs boson scales with the mixing between bino and 
higgsino components, or wino and higgsino terms) and the structure of the scalar sector is 
very constrained (for negative // however, the two contributions have negative interference 
and the scattering cross section may become small). 

From the point of view of indirect dark matter detection, the regime with bino-like 
neutralinos in the coannihilation tails is twice disfavored: first there is the S-wave sup- 
pression of the cross section for the annihilation of non-relativistic Major ana fermions into 
light fermions (at non-zero temperature, in the early Universe, the annihilation cross sec- 
tion is dominated by the P-wave) ; at the same time, the effective annihilation cross section 
((Tefff ) at freeze-out is dominated by coannihilating partner contributions, with the actual 
neutralino annihilation at T = being much smaller. The funnel regime is instead more 
interesting, since resonant annihilations would still be effective at T = 0, giving large values 
for {aesv){T = 0). Actually, most often, the enhancement effect is not at full strength, 
because the neutralino thermal relic abundance of models matching too closely the reso- 



nance condition tends to be smaller than needed: the choice falls then on models slightly 
off-resonance, again with an effect that is less dramatic for the thermally averaged term in 
the early Universe than for the annihilation cross section at zero temperature in today's 
structures. 

A similar feature emerges in the focus point region, or in MSSM models defined at the 
weak scale keeping // and Mi as free parameters, but still, say, with some common mass 
scale in the squark sector, which again is forced to be heavy by accelerator constraints. 
Requiring the thermal relic density to lie within the WMAP selected range, is equivalent, 
at a given neutralino mass, to fix the mixing between bino and higgsino components; except 
for very heavy neutralinos in the TeV range, where detection with any method looks in any 
case problematic, the mixing is large, and the annihilation rate into weak gauge bosons 
(which is not S-wave suppressed) is never at full strength, with again, to some extent, low 
temperature effects playing against indirect detection. The large mixing enhances instead 
the direct detection rate, but again this is not an effect of crossing symmetry, i.e. a feature 
which follows automatically from the definition of thermal WIMP. 

From this discussion it emerges that the hierarchy between direct and indirect detection 
may be inverted, in supergravity inspired models, relaxing some of the so-called universality 
assumptions, for instance with more complex scalar sectors [l^]. We choose here, instead, 
to restrict to minimal scenarios, and to relax, instead, the constraint coming from the 
computation of the thermal relic density. We will include in our discussion models with low 
relic densities, or, more precisely, entertain the possibility that the WIMP relic abundance is 
enhanced compared to the standard scenario. There are numerous possibilities on how such 
enhancement may take place, ranging from cosmological scenarios with faster expansion 
rates at the time of decoupling, such as in cosmologies with a quintessence energy density 



term dominating at the WIMP freeze-out temperature |19, 2C, 21], anisotropic Universes 
with effective shear density terms [p2| , |2l[] or in scalar-tensor theories ||2^ , to particle physics 
models with non-thermal sources of WIMP dark matter, such as neutralino productions 
from gravitinos or moduli decays |24 



2.2 The three benchmark scenarios 

We focus on three minimal SUSY frameworks with the lightest neutralino being respectively 
bino, wino and higgsino like, each of them in large purity configurations. All models are 
soundly motivated from the high energy physics setup point of view. The neutralino mass 
ranges within values allowed by all accelerator direct and indirect bounds. 

Binos: the Funnel model. We consider the mSUGRA scenario, and select models 
within the funnel region, along a section of the mi/2-mQ parameter space in which the 
lightest neutralino mass exactly matches the resonance condition with the CP-odd Higgs 
boson mass, i.e. 2 ■ m^ = tua- In mSUGRA, the tree level value of niA is fixed by tan /3 and 
by the sign of /i through EWSB Conditions; taking into account that we need the lightest 
neutralino to be the LSP, the resonance is realized only at tan/3 > 35 in case sgn^ < 0, 
or at tan/? > 50 if sgn/i > 0. Since for ^u < and large tan/3 the SUSY contributions to 
BR(b — > S7) are large, and rule out models up to fairly heavy neutralinos, we consider an 
example with positive fi, tan/3 = 55 and trilinear coupling ^0 = (the latter is fixed for 



Model 


Ml/2, "i3/2, ^3 


tan/3 


sgn(^) 


Defining Condition 


Funnel 
AMSB 

NUGM 


700 -:- 1450 

23 ^ 231 
879 ^ 1096 


55 
50 
50 


>0 
>0 
>0 


2 m^ ~ rriA 

mo = 1500 GeV 

M1/M3 = 10, M2/M3 = 2, H = 99.8% 



Table 1: The three SUSY models under consideration. Models Funnel and NUGM have 
mSUGRA-like boundary conditions at the high energy scale, while model AMSB features minimal 
Anomaly Mediated SUSY Breaking. A/1/2 indicates the value (in GeV) of the common gaugino 
mass at the grand unification scale Mgut in model Funnel (second column); in model AMSB the 
neutralino mass is instead driven by the value 771.3/2 of the gravitino mass (in TeV, second column); 
finally, model NUGM {non-universal gaugino masses) has the gaugino non-universality pattern 
of the 200 representation of the symmetric product of two SU{5) adjoints, and the high energy 
parameter we use is the gluino mass M3 (in GeV, second column). For each model, the value of the 
common scalar SUSY breaking mass mQ is dictated by the Defining Condition, which is indicated 
in the last column of the Table. The scalar trilinear coupling Aq ^ for all models. 



Model 


Bino fraction 


Wino fraction 


Higgsino fraction 


Funnel 
AMSB 

NUGM 


> 99.6% 

< 0.02% 

< 0.01% 


< 0.05% 

> 98% for m^ > 100 GeV 

< 0.2% 


< 0.5% 

< 2% for m^ > 100 GeV 

99.8% 



Table 2: The lightest neutralino composition for the three SUSY models of Tab |l|, in terms of the 
bino, wino and higgsino fractions. 



definiteness, as it does not play much of a role). The resonance curve in this configuration 
starts at a minimum value of mi/2 (and therefore ttt,^) below which the lightest stau is 
lighter than the lightest neutralino^, and ends at at maximum mi/2 corresponding to the 
model with neutralino relic abundance exceeding the CDM upper limit. The bino purity 
of these models is anywhere significantly large, always above 99.6% (see tab. ^). We will 
hereafter indicate this set of configurations as the Funnel model. 

Winos: the AMSB model. In the so-called minimal Anomaly Mediated SUSY Break- 



ing (mAMSB) scenario |25, 26, 27, 28] the gaugino masses are proportional, through the 
gravitino mass 77^3/2, to the beta functions of the SU{3), SU{2) and ^(1) gauge groups, 
therefore drastically affecting the soft breaking gaugino mass hierarchies at low energies, as 
compared to a supergravity and grand-unification inspired scenario: the lightest neutralino 
is in fact always wino-like, with a remarkably large purity. This feature induces a nice prop- 
erty of the mAMSB scenario, i.e. the homogeneity of observables over its parameter space: 
except for accidental cancellations, the relic abundance and the rates for indirect detection 
are sensitive to the value of the LSP mass, but have a rather mild dependence on other 
specific features of the low energy structure of the theory. We decided to resort to a rather 
large value of the common scalar mass parameter mg = 1500 GeV, and of tan/3 = 50. The 
sign of ^ has been assumed to be positive. The range of the gravitino mass (and therefore 
of the lightest neutralino) is bounded from below by the null-results of chargino searches at 



^In this low mass region, stau coannihilation effects are also present, and are taken into account in the 
computation of the relic abundance; they are, however, not relevant in the context of dark matter searches. 



LEP [29p, and from above by the fact that one enters in the region where EWSB can no 
longer be fulfined, giving raise to an unphysical Higgs sector. We label the configuration 
defined by this choice of parameters, which is summarized in tab. ||, as AMSB model; 
details on the composition of the lightest neutralino are given in tab. ^. 

Higgsinos: the NUGM model. The assumption that gaugino masses unify at the 
GUT scale Mqut niay be relaxed in a number of ways. In the context of supergravity, non- 
vanishing gaugino masses are generated through the SUSY breaking vev of the auxiliary 
component of the lowest order non-renormalizable term in the gauge kinetic function. The 
SUSY breaking vevs, in order to preserve gauge invariance, must lie in a representation 
belonging to the symmetric product of two adjoints of the underlying unified gauge group. 
In the case of SU{5) GUTs, one has 

(24 X 24)^y^„^ = 1 e 24 e 75 e 200. (2.2) 

Universal gaugino masses are generated only if the SUSY breaking field lies in the singlet 
representation 1. In all other cases, particular ratios between gaugino masses will hold at 
the GUT scale [30, ^]. In particular, two of the four representations yield a higgsino-like 



LSP, namely the 75 and the 200. We focus here on the 200, because in the case of the 75 
successful EWSB forces the parameter space to rather narrow regions, and to low values 



of tan/3 |32|. The parameter space in the corresponding minimal scenario is analogous to 
the mSUGRA case, with a gaugino mass parameter, which we take to be the gluino mass 
M3, and a common scalar mass parameter ttt-q (plus again tan/3, Aq and sgn(/x)). Since 
we want to study the case of a pure higgsino, we solve the parameter space degeneracy 
by requiring a fixed higgsino content of the lightest neutralino, H = 99.8%. This specific 
value has been chosen in order to maximize the neutralino mass range, rriy. turns out to 
decrease at growing mo, along the iso-higgsino content lines; analogously to our previous 
choices, we fix tan/3 = 50, Aq = and sgn/i > 0, and let ttiq vary from to the largest 
value compatible with EWSB, achieving a rather wide neutralino mass range. Details on 
the model are provided in the bottom line of tab. |l|; this configuration will be indicated 
from now on as NUGM (non-universal gaugino mass) model. 

For all the three benchmarks scenarios introduced above, soft breaking parameters, 
gauge and Yukawa couplings are evolved down to the weak scale with the ISASUGRA 
RGE code as given in version 7.67 of the ISAJET software package [^] (introducing some 
minor changes, such as, the conversion of the ISASUGRA code to double precision to 
improve on its stability; for more details on this and other technical points on the code im- 
plementation and on the interface with the DarkSUSY package, see the discussion in \p%). 
Weak scale spectra are interfaced into the DarkSUSY computer code [p5|| , which is then 
used for computations of the relic density and of direct and indirect detection rates. Re- 
garding, in particular, the relic abundance calculations, we rely on a code which provides 
a high-precision fully-numerical treatment of resonances and coannihilation effects, here 



^In the case of mAMSB, the LEP bound on the chargino mass is somewhat weakened due to the very 
small splitting between the lightest neutralino and the lightest chargino. We will nonetheless restrict to 
values larger than 100 GeV. 
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Figure 1: In the figure to tiie left we sliow tlic neutralino relic abundance, in the standard cos- 
mological scenario, for the three models described in the text (see Tab. fl), respectively featuring 
a resonantly annihilating bino (Funnel), a wino (AMSB) and a higgsino (NUGM) lightest neu- 
tralino, as a function of the neutralino mass. The green band indicates the preferred WMAP range 
|3|] . The right panel shows instead the relic density enhancement within the quintessential scenario 
of ref. ||36| , i.e. with an exponential potential, with the exponent A = 3.5, for the three benchmark 
models. The upper part shows the relic abundance, while the lower part the relative enhancement 
AJ7 — (J7q — r2)/r2, both as a function of the parameter ^0, the ratio of the quintessential and the 
radiation energy densities at the neutralino freeze-out temperature without quintessence. 



necessarily needed to properly include the resonance on the A Higgs boson for the Funnel 
model and chargino coannihilations for the AMSB and NUGM models (for both pure 
winos and higgsinos, the lightest chargino is nearly degenerate in mass with the lightest 
neutralino). Neutralino relic abundances for the three models considered, for a standard 
cosmological setup and no extra non-thermal sources, are shown in the left-hand side of 
Fig. y, in the neutralino mass range allowed for each scenario. 

In the right panel of Fig. || we show instead an example of the quintessential enhance- 
ment of the thermal neutralino relic abundance at work on three models within the above 
described benchmark scenarios, at a common mass m^ = 300 GeV. The upper panel shows 
the relic abundance of each model as a function of the parameter ^0, which is defined 
as the ratio between the quintessential energy density p^ (most likely, in a "kination " 
phase, i.e. a phase of fast-rolling and p^f, scaling like a~^, with a the scale factor of the 
Universe) and the radiation energy density /Jrad at the neutralino freeze-out temperature 
in the absence of quintessence. This parameter has been shown to suitably gauge out most 
of the dependence on the parameters in the particle physics setup, and to describe the 
relic density enhancement effect in other cosmological scenarios, as for instance the case 
of an anisotropic primordial Universe with an effective shear energy density scaling again 
as a~^ |22, 21|. From the figure we deduce that, at a neutralino mass of 300 GeV, the 



Funnel model would require an extra energy density component 10 times larger than that 
of radiation at neutralino freeze-out, while larger amounts are needed in the NUGM and 
in the AMSB cases (respectively around 10^ and 10'^). The lower panel shows the relative 
enhancement, 

An = -^ — . (2.3) 

As already mentioned and described in further details in ref. |20|, we find a rather small 
spread in values of AQ among the SUSY models considered here and for given values of 
the parameter ^^. 

3. Indirect searches through antimatter cosmic ray measurements 

Since the first proposals to search for the exotic antimatter components in cosmic rays due 
to dark matter WIMP pair annihilations [ p7| , [3^ ], calculations of the expected fluxes have 
been performed with an increasing degree of sophistication. In recent estimates, see, e.g., 
refs. |39, HH, 41, 42, 43|, refinements have regarded both the modeling of source functions 



and the description of the propagation of charged cosmic rays in the Galaxy. 

3.1 Estimates of the antimatter fluxes 

We discuss here the case of neutralino-induced antiproton, positron and antideuteron cos- 
mic ray fluxes. Predictions for the three species are derived in a consistent framework after 
going through a few steps. As already stressed, source functions are proportional to the 
annihilation probability for a pair of WIMPs, i.e. {aef[v){T = 0), and depend on neutralino 
mass and on the branching ratios for the various annihilation channels. Annihilation yields 
fragment and/or decay, and generate stable antimatter species; for positrons and antipro- 
tons, the simulation of this process is performed using the Pythia [44| 6.154 Monte Carlo 



code, in each of the allowed two-body final state and for a set of 18 neutralino masses. For 
D sources we implement instead the prescription suggested in ref. |42] to convert from the 
p-n yields. 

Source functions are proportional as well to the number of WIMP pairs which poten- 
tially may annihilate at each place in the Galaxy; postulating that dark matter is smoothly 
distributed in the Galactic halo, the number density of pairs can be expressed in terms of 
the dark matter density profile as 1/2 (p^(x )/r7i^)^. Even in this limit, further extrapo- 
lations are needed, as, unfortunately, Px{x) is poorly constrained by available dynamical 
data on the Galaxy. In this analysis, we will mainly focus on a dark matter halo described 
by the so-called Burkert profile ||45[|: 



(1 + r/a) (1 + {r/aY) 
it is a profile with a large core radius, that has been tested against a large sample of the 



rotation curves for spiral galaxies |46]. Such choice is very conservative, as results from 
N-body simulation of hierarchical structures in CDM cosmologies find instead singular halo 
profiles. In an extreme model for the baryon infall [E^], with very large angular momentum 



transfer between the baryonic and the dark components, a Burkert-type profile may replace 
the CDM cuspy halo after the gas has been settling in the inner portion of a galaxy to form 
its luminous components. Scenarios of baryon infall in which the CDM cusp is preserved 
are feasible as well, and we will discuss one such possibility in sec. |^, mentioning also the 
implications for neutralino dark matter detection. We make a sample choice of the free 



parameters in eq. 3.1, fixing the length scale parameter a = 11.7 kpc and the local halo 
density pB{ro) = 0.34 GeV cm~^; this configuration has been found after implementing 
available dynamical constraints, including constraints on the local halo density from the 
motion of stars in the Sun's neighborhood, on the total mass as derived from the orbit 
of the outer satellites, as well as by requiring consistency with the Milky Way rotation 
curve and with measures of the optical depth towards the galactic bulge. We have also 
taken into account the correlation between halo mass and halo concentration as found in 
the numerical simulations (in the structure formation language, the halo we consider has 
a virial mass and concentration parameter of, respectively, Mmr = 1-3 x lO^^M© and 
Cmr = 16, see ref. [48| for details). 



Once the source functions are fully specified, the next step is to model the propaga- 
tion of charged cosmic rays through the Galactic magnetic fields. We consider an effective 
two-dimensional diffusion model in the steady state approximation. We do not include ex- 
plicitly reacceleration effects in this propagation model, but mimic them through a diffusion 
coefficient which takes the form of a broken power law in rigidity, R, 

D = Do {R/Rof-^ if R>Ro 

D = Do a R<Ro. (3.2) 

This form has been used in a number of studies on the propagation of cosmic rays; e.g., in 
ref. 1^^, using the Galprop [|50|| propagation code, it has been shown that it gives a fair 
estimate of ratios of primary to secondary cosmic ray nuclei, for a suitable choice of the 
free parameters in the model. We implement here this same setup, i.e. take eq. ( |3.2D with 
Dq = 2.5 X 10^^ cm^ s~^ and Rq = 4: GV, in a cylindrical diffusion region of radius equal 
to 30 kpc and half height equal to 4 kpc, plus a galactic wind term. For antiprotons and 
antideuterons, this setup is interfaced to the semi-analytic diffusive-convective propagation 



model described in ref. [40|, which do not allow for an energy loss term (particles are 
removed whenever they scatter on the interstellar medium). For positrons, we exploit the 
results in ref. [^], where a propagation model with a term accounting for positron energy 
losses (the dominant terms accounting for inverse Compton scattering on starlight and 
the cosmic microwave background) was solved analytically; such model has been recently 
improved and extended to allow for the implementation of a diffusion coefficient in the 
form of eq. ( |3.2[ ), and to keep a full two-dimensional structure (the corresponding code is 
included in the latest release of the DarkSUSY package). 

The final step to compare with measurements is to include a treatment of solar modu- 
lation effects, i.e. the propagation of cosmic rays through the solar system up to the Earth 
location, against the solar wind. To sketch this effect, we implement the one parameter 
model based on the analytical force-field approximation by Gleeson &; Axford |51] for a 
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Figure 2: The interstellar antiprotons flux, as a function of the antiprotons kinetic energy Tpbar, 
for two values of the neutralino mass, respectively m^ = 300 GeV (left) and m^ = 500 GeV (right), 
for the three SUSY models of Tab. 0. 



spherically symmetric model. The solar modulation parameter, sometimes dubbed Fisk 



parameter ^p |52|, is for simplicity assumed to be charge-sign independent, and it is as- 
sumed that, for all species, it can be chosen to be equal to the corresponding parameter as 
determined for proton cosmic-ray flux. The latter has been measured with some accuracy 
over the full 11-year solar cycle. 

3.2 Antimatter spectral features for the three benchmark models 

Antimatter fluxes from neutralino annihilations depend on one hand on the number density 
of neutralinos in the galactic halo (and therefore, for a given CDM halo profile, on the 
inverse of the neutralino mass squared), and on the other hand on the specific features 
of the annihilation processes of neutralinos, in terms of the overall interaction rate and of 
the standard model final products. The three benchmark models considered in the present 
work feature relatively large annihilation cross sections, when compared to the rule of 
thumb estimate, approximately valid in a standard cosmological scenario, 



(<7eflff) 



3 • 10^27 ^^3 g- 



(3.3) 



The values of {aegv), at a sample neutralino mass of 300 GeV, lie in fact between 10~^^ 
and 10~^^ cm^s~^. As expected, low relic density models will thus have noticeably large 
antimatter fluxes, as compared to models which give the expected relic abundance within 
the standard cosmological scenario. 

We show in figures y, y and Q the interstellar fluxes of antiprotons, positrons and 
antideuterons for a few selected models within the three benchmark scenarios. In the left 
and right panels of each of the figures we compare (kinetic) energy spectra for the given 
antimatter species at different neutralino masses, respectively 300 GeV (left) and 500 GeV 
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Figure 3: The interstellar positrons flux, as a function of the positrons kinetic energy Tepius, for 
two values of the neutralino mass, respectively m^ = 300 GeV (left) and m^ = 500 GeV (right), 
for the three SUSY models of Tab. 0. 

(right). The l/rn^ suppression effect due to neutralino number density dilution at growing 
masses is everywhere apparent. 

The spectral features of the antimatter fluxes, on the other hand, critically depend 
on the decay channels of the annihilating neutralinos. In the case of the Funnel model, 
the largely dominating decay channels are 66 (around 90%) and t~^t~ (around 10%): this 
pattern is due to the fact that the value of tan /3 which has been chosen is particularly large, 
so that decays into up-type fermions (charm or top) are very suppressed with respect to 
those into down- types. The resulting antimatter from hadronization of 6 quarks as well as 
r decay products yield typically soft antimatter particles: as emerging from the figures, at 
the lowest energies of interest, binos in the Funnel model always give larger fluxes with 
respect to higgsinos and winos. 

The energy spectrum in the AMSB model reflects the large degree of purity of the 
dominating decay channel, namely that in W~^W~ , whose branching ratio is always above 
98%. As regards hadronic antimatter (antiprotons and antideuterons) , annihilating winos 
tend to give rise, from hadronization of W decay products, to a peak in the antimatter 
particle energy around a few GeV. The precise location of the peak non-trivially depends 
on the details of the antiproton and antineutron production, though its location in energy 
increases at larger masses (see fig. |2| and ^. In the positron channel, a clean peak is located 
at Sg+ ~ m^/2: the origin of this very neat feature is the prompt decay W^ — > e~^ fg, 
where, since the annihilating neutralinos are non-relativistic, the energy of the decaying 
W is close to m^. 

Turning to the case of higgsinos (NUGM model), the spectral features of antiprotons 
and positrons neatly reproduce those of winos, though at lower fluxes (the suppression being 
motivated both by a smaller effective annihilation cross section and by a suppression factor 



12 



10 



■S 10 - 



E 
I 10"' 



: 10" - 



10" 



Funnel 

AMSB 

NUGM 



m =300 GeV 



- ■■.. \ 
^ ■•., \ 



10 10 

To»„ [GeV] 



\ 

,'■. \ 

\\ \ 



\\ \ 
\ \ 



10' 



^10"' 

o 

0- 
§10-' 



m =500 GeV 



10"" - 



10- 



Funnel 

AMSB 

NUGM 



10 10 

To.„[GeV] 



\ 

■. \ 
-■■.. \ 
v.. \ 



\ 



\ 

. \ -I 

\\ 

\ 



10' 



Figure 4: The interstellar antideuterons flux, as a function of the antideuterons kinetic energy 



Tbbar, for two valucs of the neutralino mass, respectively m^ 
(right), for the three SUSY models of Tab. |l|. 



300 GeV (left) and m^ = 500 GeV 



in the couplings of higgsinos to the W bosons) . Noticeably, the antideuteron flux is instead 
larger than that of winos at energies below the GeV. This pattern is motivated by the fact 
that higgsinos decay with ~ 50% probability into W~^W~, and with ~ 45% probability 
into Z^Z^: while antiprotons and positrons production from Z^ decays is very similar to 
that from W decays, the creation of low energy antineutrons from Z^Z^ is remarkably 
enhanced, thus giving rise to the mentioned peak in the antideuteron low-energy flux. 
Fig. S shows the spectral features, after solar modulation for a given step along the 



solar activity cycle (see the discussion in sec. 3.1), of primary antiparticles and of the 
background, comparing the total expected signals to the data on antiprotons and positrons 
taken during the corresponding modulation phase. The figure refers to a common neu- 
tralino mass of 300 GeV. As regards the secondary antimatter fluxes, which play here the 



role of backgrounds, our estimates are produced running the Galprop |50| code in the con- 
figuration for propagation parameters we have adopted for the signals. We remark that for 
both species, the computed backgrounds provide by themselves excellent fits of the data: 
we obtain, for background only, a reduced x^ equal to 0.82 for antiprotons and to 0.95 for 
positrons. A primary flux will then be excluded as long as it gives rise to a statistically 
unacceptable x^i iio longer compatible with currently available data. At a given neutralino 
mass, the signal-to-background ratio for antiprotons, after solar modulation, is much more 
promising than that of positrons, with the possible exception of large energies, around the 
m^/2 peak in the positron spectrum, which lies however far above the currently sampled 
energies. 



4. Current exclusion limits 

Ruling out a model from antimatter fluxes amounts to evaluate the x^ of the expected 
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Figure 5: (a): The solar modulated antiprotons flux, as a function of the antiprotons kinetic 
energy Tpbar- The black line corresponds to the calculated background, while the three colored 
thick lines to the total signal for the three SUSY models at a mass ra^ = 300 GeV. The thin lines 
correspond to the SUSY contributions alone. The data from BESS-98 || and CAPRICE-98 111 
are also shown, (b) The solar modulated positrons flux, as a function of the positrons kinetic energy 
Topius- The black line corresponds to the calculated background, while the three colored thick lines 
to the total signal for the three SUSY models at a mass m^ = 300 GeV. The thin lines correspond 
to the SUSY contributions alone. We also plot the positron data from MASS-91 ||3l, HEAT-94/95 
155|1 and CAPRICE-98 fs^. 



signal plus background, after taking care of the different solar modulation effects in the 
relevant period of data-taking. We will show here the reduced x^, i-e. the x^ divided by the 
relevant number of data (respectively 49 for antiprotons and 32 for positrons). The 95% 
confidence level (C.L.) exclusion limit will lie at a x^ around 1.33 for antiprotons and 1.4 
for positrons, though x^ ~ 1 indicates that the model is around its ^^visibility threshold!\ 
We will find models for which reduced x^ ^^^ above 10 or so, which are excluded with no 
doubt, and models closer to the quoted confidence levels; for the latter one should keep 
in mind that we are not taking into account uncertainties in the propagation parameter 
model and in the halo profile (although in this respect we are taking a rather conservative 
scenario), hence limits we will quote should not be intended as sharp cut-offs. 

The quality of the data on the local antiproton and positron cosmic ray flux has kept 
improving in recent years. As regards the antimatter flux data we take into account, we 
will compare the predicted fluxes to the antiproton data collected by the BESS experiment 



during its flights in 1997, 1998, 1999 and 2000 [§3|, with fairly good statistics, in the 
energy range between 180 MeV and 4.2 GeV, and by the CAPRICE experiment during its 
1998 flight [^4[ in the range between 3 and 50 GeV. For the positron fluxes, we consider 



the data pubhshed by the HEAT Collaboration about the 1994-1995 flight |5|], by the 
CAPRICE team in a flight in 1994 ^, and finally data obtained by MASS-91 in 1991 
1 57]: the overall energy range covered by these measurements extends from 460 MeV to 
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Figure 6: Current exclusion limits on the AMSB model from direct and indirect Dark Matter 
searches, as a function of the neutralino mass. In the upper panel we show the reduced y^ of 
the background plus SUSY signal, and, as a guideline, x^^^ = 1; in the lower panel we show the 
ratio of the expected signal and of the current exclusion limits ( Visibility Ratio) for direct dark 
matter detection (neutralino-proton spin-independent scattering cross-section, solid line) and on 
the neutrino-induced muon flux from neutralino pair-annihilation in the center of the Sun (dotted 
line). The ^^ Visibility Line" corresponds to values of the expected signal equal to the current 
exclusion limits: models above the Visibility Line would be already ruled out. 

34.5 GeV. We have chosen not to include in our analysis data which have been reported 
just as antiproton or positron fractions (rather than absolute fluxes) and datasets such as 
the one on positrons from the AMS test flight |5^] mapping a low energy interval in which 
a primary neutralino-induced contribution is expected to be rather suppressed. 

We compare the statistical analysis on antimatter fluxes with the expected signal- 
to-sensitivity {^^visibility") ratio for direct, spin-independent searches and for the muon 
flux at neutrino telescopes originating from neutralino annihilations in the center of the 
Sun'^. The visibility ratio for direct detection refers to the expected neutralino-proton 
scattering cross section over the current experimental sensitivity level. For definiteness, we 
will refer to the EDELWEISS 2002 [^] results and estimate the corresponding exclusion 
limit, rather than within the standard Maxwell-Boltzmann local velocity distribution as 
usually done, by using the self- consistently derived velocity distribution of the halo model 



^We always find that spin-dependent direct detection, and the muon flux from the center of the Earth 
or from the Galactic center have visibility ratios by far smaUer than those we consider. 
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Figure 7: The same as in fig. ^ again in the AMSB model, but with a rescaled neutralino density 
profile, according to the formula p'^'^'^ = pcnm ^coM/Min (Jl^ 



■■X' ^^cdmJ 



under consideration, and taking into account relevant effects, such as target materials, form 
factors, and threshold. ^ 

As far as the muon flux induced by neutrinos from the sun, the current best limits 
are from the SUPER-KAMIOKANDE Collaboration in 2002 ^. Also for this detection 



method, the signal is computed estimating capture rates with the appropriate velocity 
distribution. 

4.1 The AMSB model 

Fig. ^ shows the current visibility ratios and exclusion limits for the AMSB model. Re- 
markably, we find that direct detection, in this particular case, is the less promising search 
strategy, with expected signals lying more than two orders of magnitudes below current sen- 
sitivity. Neutrino telescopes do little better. On the other hand, masses below ~ 270 GeV 
are ruled out, within the conservative Burkert halo model we consider, from the overall fit 
to antiprotons data, while the expected positron flux is above current sensitivity for masses 
below approximately 170 GeV. In case we had assumed that the model under scrutiny is 
describing some suhdominant component of the full CDM term, and if we had accordingly 



*An effect compatible with being due to WIMP-nucleon scatterings has been reported by the DAMA 
Collaboration, see their latest report in ref. fed] . Unfortunately, none of the models considered here give an 
effect at the level of such signal. 
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Figure 8: Direct spin independent neutralino-proton scattering cross section (left) and Muon flux 
induced by neutrinos produced by neutralino annihilations in the center of the Sun (right) for 
various (around 2500) minimal Anomaly Mediated SUSY breaking models. We linearly scanned 
2 < tan/3 < 50, 200 < mo < 5000 and lO'^ < TO3/2 < 4- 10^ and picked with equal probability either 
sign of fi. Four particular parameter space slices at rriQ = 1500 GeV, and at different tan/3 and sign 
of /i are also shown as guidelines. 

rescaled the neutralino densities according to the thermal rehc abundance computed in a 
standard cosmological scenario, i.e. applying the rescaling recipe: 



nr 



PCOM 



Min([7^,[7Si,°^) 



(4.1) 



results would be largely altered, as we show in fig. |^: direct detection and neutrino tele- 
scope visibility ratios would have been both five orders of magnitude below one, and an- 
tiprotons and positron fluxes would give rise to a signal completely indistinguishible from 
the expected background. Relic density enhancement mechanisms are therefore mandatory 
within the minimal AMSB framework to achieve signals compatible with current, as well 
as with future experimental sensitivity. We find that the same conclusion applies for the 
other benchmark models we consider here. 

As manifest from the shape of the corresponding curve in fig. ^, direct detection rates 
are particularly low because a cancellation in the contributions from the i-channel ex- 
changes of Hf and H2 takes place in the particular parameter slice we consider here as 
a benchmark scenario. For clarity, we performed a random scan of the full mAMSB pa- 
rameter space (see the caption in fig. ^ for details on the scan). We report in the scatter 
plots of fig. ^ our results, together with three parameter space slices as guidelines and 
with the considered AMSB model as well. The left part of the figure refers to direct 
detection, while on the right we display muon fiuxes at neutrino telescopes. Notice that, 
although larger scattering cross sections may be obtained, all models are far below current 
sensitivity for direct detection, therefore not altering our conclusions. As regards neutrino 
telescopes, though very few points lie above the current exclusion limits, we point out the 
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Figure 9: Current exclusions limits on the Funnel model. Symbols and conventions are the same 
as in fig. O. 



the resulting masses are always smaller than the limits we derived from current data on 
antimatter searches. 



4.2 The Funnel and the NUGM models 

Current exclusion limits for the Funnel and for the NUGM models are showed in figures 
pi and 10 respectively. In the case of the Funnel, where masses are rather large, we do not 
find any point which is neither ruled out by current data on antimatter fluxes nor by direct 
detection and neutrino telescopes. The latter two experimental techniques have achieved 
so far a sensitivity three orders of magnitude or worse below the expected signals. 

Regarding the NUGM model, once again we single out a case for which direct dark 
matter detection is the least competitive detection strategy. Neutrino telescopes, though 
featuring a one order of magnitude better visibility ratios, are also a factor 100 or more 
below the needed sensitivity. Turning to antimatter, in the NUGM scenario positrons do 
not give any statistical constraint, while antiproton expected fluxes are not compatible with 
current data for low masses. We remark that this is another example of a SUSY model 
which is currently constrained only by antimatter searches; moreover notice that this is 
true with the rather conservative Burkert halo profile considered here, and that constraints 
from antimatter searches, within other profiles, can be much tighter (see sec. ^). 
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Figure 10: Current exclusions limits on the NUGM model. Symbols and conventions are the 
same as in fig. |[ 

5. Prospects for SUSY DM detection through antimatter searches 

5.1 Discrimination perspectives at future antimatter search experiments 

New generation space-based experiments for antimatter searches PAMELA [62| and AMS [63] 
will tremendously enhance the resolution and accuracy of positron and antiproton spectra 
measurements, as compared to existing balloon borne results (see sec. y). With the purpose 
of assessing discrimination capabilities of future experimental facilities, we will sketch here 
the possibility of disentangling an exotic component out of a standard secondary back- 
ground. To this extent, we will implement a statistical x^ analysis to compare the case of 
a pure background measurement to that of the occurrence of a signal. 

The relevant experimental parameters entering the estimate are given by: 

- The geometrical factor of the experimental facility, i.e. its effective area. A] 

- The time of data acquisition T; 

- The energy coverage of experiment, with the relative definition of energy bins, i.e. 
its number nb and the size {AE)i of each of them. 

We will declare that a given SUSY model is going to be discriminable at a certain 
future experiment, and at given X% confidence level, if the x^ induced by the SUSY model 
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is larger than the (x^)m, corresponding to nb degrees of freedom. Letting Nf = N^ +N^ 
be the number of projected events in a given bin i, the sum of the number of signal plus the 
background events, and iV^P be the number of observed events, with a standard deviation 
A »fO , the x^ is defined as 

i 

X' = f: ^''] ~ %^\ Nr = N^ + Np. (5.1) 

We will suppose that the standard deviation has a Gaussian distribution, i.e. 

^No ^ 4nO. (5.2) 

We are interested in finding the limiting cases, i.e. those cases for which an eventual 
signal is a small component with respect to the background. We will therefore make the 
assumption that: 

N^ « iVf , or A^P :^ iVf . (5.3) 



Eq. (|5.1| ) will then read 




Now, since the number of events in an energy bin /S.E is given, as a function of the flux of 
particles (/>, by 

N = {AE) -(P-A-T, (5.5) 

and indicating with (ps and (ph the signal and background fluxes respectively, eq. (|5.4| ) will 
read 

X' = E ^■{AE),.A.T. (5.6) 

The quantity in eq. ^^ is what will be used to asses the future sensitivity at antimatter 
experiments. We will declare that a model is within discrimination capabilities of a given 
future experiment at X% confldence level if it satisfies the relation 

X' > (X^)^7". (5.7) 

We focus, for definiteness, on the case of the PAMELA detector, and compute the 
reduced x^ for an effective area of 24.5 cm^sr, an exposure time of 3 years, and resorting to 
a trial energy binning as sketched in ref. 163] ^. The results we will show are in the limit of 
known background, i.e. in the (optimistic) scenario in which degeneracies in the parameters 
used to model the propagation of charged cosmic rays in the Galaxy are resolved, say, by 
precision measurements of ratios of secondaries to primaries for several light cosmic-ray 



^The data binning is not going to be homogeneous, since the fluxes have a peak at a few GeV energy. 
For definiteness, we divided the total logarithmic energy interval scanned by Pamela into three parts and 



estimated the respective number of bins from fig. 7 of ref. [p4|. 



20 



nuclei; in this context, we model the background according to the same estimates already 
implemented in the previous section. 



Fig. 11 shows the future discrimination perspectives for the three benchmark models 



described in sec. 2.2 for antiprotons (left) and positrons (right). The horizontal green line 
represents the 95% C.L. exclusion limit for PAMELA, while the three lines indicate the 
quantity in the left-hand side of eq. ( |5.6| ), i.e. the projected x^ for the given SUSY model. 
Notice that the two horizontal lines are slightly different, as the number of bins nb for 
antiprotons is larger than that for positrons. 

For wino and higgsino-like lightest neutralinos, antiprotons will probe a larger portion 
of parameter space than in the case of positrons. On the other hand, a bino-like neutralino 
will give a signal which is not going to be detected even in the favorable case of the 
occurrence of an s-channel resonance. 

We wish to make another step forward and present the results we just obtained in 
a form which is less closely tied to what it has been assumed as properties of a given 
detector. In fact, we definitely do not aim to assess in detail the detector responses, and 
what we extrapolated for PAMELA should be taken just as an approximation which can 
be improved with dedicated simulations. We wish instead to present a simple scheme to 
which an eventual refined analysis can refer to, without having to readdress from scratch 
the problem for the particle physics point of view; at the same time, this should be a 
valuable tool to make it easy to extend our analysis to other alternative particle physics 
scenario. With this goal in mind, we notice that the upcoming antimatter experiments 
PAMELA and AMS will collect a rather large number of positrons and antiprotons, and 



hence the number of energy bins in eq. |5.6| will be significantly large. We find therefore 
appropriate to approximate the discrete sum in this equation with an integral. The region 
of integration will be given by the lowest and the largest energies, Tmin and T^aax, accessible 
by the experiment. We henceforth define the quantity 



The visibility condition, which reflects the continuum version of eq. |5.6| , is then given by 

I^ > ^^- (5.9) 

In the parameter J^, the dependence on the extreme of integration is small as long as 
the peak on the signal to background ratio falls within them (and this is always the case 
for the models we are considering and for values expected for PAMELA and AMS); we 
have assumed that the effective area is independent of energy, while, if this is not a good 
approximation, a weight function should be accordingly introduced. For the PAMELA 
setup already exploited, supposing a number of bins given by n^ = 55 -^ 60, and taking 
(x^)nu ~ "^5, we find the following critical values for J^, after, respectively, one and three 
years of data taking: 

/]y = 9.7 . IQ-s cm"2g^-ig-i /3y = 3.2.10"'^ cm"2gj.-ig-i^ (510) 
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Figure 11: Discrimination perspectives on the three SUSY models of Tab. n^ for the PAMELA 
experiment, as regards antiprotons (left) and positrons (right). The three lines correspond to the 
expected total x^ for the three SUSY models, calculated assuming a small signal to background 
ratio (i.e. in the approximation of eq. ( |5.3| )). The horizontal solid line reproduces the 95% C.L. 
exclusion limit Xq 95 foi' t^^ assumed PAMELA data binning. 



We will refer to these values when comparing the perspectives of detection with different 
methods. 

An ambitious goal for the future is to reach the experimental sensitivity to detect 
the antideuteron component in cosmic rays. The procedure we will follow to assess the 
detection perspectives of neutralino dark matter in this channel will be different from what 
we have just done for antiprotons and positrons. For antideuterons the background is 



expected to be totally negligible in the low energy regime [43|, and even detection of 1 
low-energy event would imply discovery of an exotic component. We refer, as the ultimate 
reach experiment in the future, to the gaseous antiparticle spectrometer (GAPS) [^. This 
is a proposal for an instrument looking for antideuterons in the energy interval 0.1-0.4 GeV 
per nucleon, with estimated sensitivity level of 2.6 x 10~^m~"^sr~^GeV~ s~^, to be placed 
either on a satellite orbiting around the earth or on a probe to be sent into deep space. 



Fig. 12 shows the visibility ratios (i.e., the projected signal over projected sensitivity) for 
GAPS-like experiments placed on a satellite around Earth (left) and on a deep space probe 
(right). 

Remarkably, antideuterons searches in the low energy range will very efficiently probe 
the funnel region: the expected signal for the Funnel benchmark model is in fact in 
this case more than one order of magnitude larger than that from winos or higgsinos. 
This enhancement of antideuteron production can be understood looking back at fig. H: 
antideuterons are efficiently produced by resonantly annihilating binos, which at large tan /? 
mainly decay into bb. Soft jets from b decay then induce, by hadronization, the low energy 
tail of antiprotons and of antideuterons (see also fig. ||). In the case of winos and higgsinos, 
instead, the main contribution comes from W decays, and features a peak at larger energies. 
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Figure 12: Discrimination perspectives on the three SUSY models of Tab. H for a GAPS-hke 
experiment for antideuteron searches in the antideuteron kinetic energy interval 0.1 GeV < T-g- <0.4 
GeV, respectively for a satellite around Earth (right) and for a probe in deep space (left). On the 
y axis we plot the expected signal-to-sensitivity ratio; the solid lines correspond to a signal from 
SUSY contributions which equals the assumed experimental sensitivity. 



ranging from 1 to 10 GeV. At large enough neutralino masses, higgsino like neutralinos 
tend to produce soft antideuterons from Z^ decays, as already pointed out from the left 
panel of fig. 0, a feature which explains the anomalous rise in the visibility ratio for the 
NUGM model in fig. 0. 

5.2 Future search strategies and the role of antimatter 

In order to compare the sensitivity of future direct dark matter searches and searches with 
neutrino telescopes with that of antimatter, we will hereafter deal with the quantity J^, 
rather than with the x^i &i^d plot the experimental discrimination sensitivity of PAMELA 
at different data taking times (see eq. ( 5.10| )). As regards, instead, antideuterons, we will 
keep on with the visibility ratio; we plot it together with the visibility ratio for direct 
spin-independent searches, assuming as reference sensitivity for future experiments that of 
the proposed XENON detector |^^, and for the muon-induced flux from the center of the 
Sun at neutrino telescopes, implementing the projected sensitivity of the km^-size detector 
which is being built by the ICECUBE Collaboration [gTJ. 



Fig. 13 compares future detection perspectives in the AMSB benchmark model. Re- 
markably, the most promising detection strategies, as emerging from the upper panel of 
the figure, reside in antiproton searches, which in one year of data-taking will probe wino 
masses up to approximately 370 GeV. For comparison, direct detection at XENON will be 
able to detect masses not larger than 200 GeV, while this model will not give any signal at 
Neutrino Telescopes. Let us stress that what we find holds quite independently from the 
value of rriQ. This in turns implies that for large values of mo, where visibility at the LHC is 
going to be much suppressed, owing to a very heavy squark spectrum, antimatter searches 
will be the only way to probe this kind of anomaly mediated SUSY breaking models. We 
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Figure 13: Future exclusions limits on the AMSB model from direct and indirect dark matter 
searches. In the upper panel we plot the discrimination parameter /^ respectively for antipro- 
tons (long dashed line) and positrons (dot-dashed), as well as the PAMELA experiment projected 
sensitivity after 1 and 3 years of data taking. The lower panel displays Visibility Ratios (expected 
signal over projected future exclusion limits) respectively for direct neutralino searches {a^-p, spin- 
independent) at XENON- 1 ton ||66| , for indirect neutralino detection through the neutrino- induced 
muon flux detection from the center of the Sun at ICECUBE [p^, and for antideuteron searches 
with a satellite around the Earth. 



point out again that this conclusion holds in case these dark matter candidates are pro- 
viding the bulk of the CDM in the galaxy (i.e. in presence of some mechanism to enhance 
the relic density), otherwise, rescaling fluxes, the discrimination capability is washed out. 
As regards the Funnel benchmark model, see fig. ^ antideuterons searches give the 
lion's share in detection perspectives, as previously pointed out. For this model, both 
antiprotons and positrons will not be able to probe much of the parameter space, since 
neutralino masses lie mostly above 300 GeV, but searches for low energy antideuterons 
might be able to reach significantly heavier masses, up to 450 GeV. A comparison with 
standard direct detection strategies, whose discrimination capability is going to be less than 
350 GeV, shows once again the utmost importance of antimatter searches in the quest for 
SUSY dark matter. 



Finally, fig. 15 reproduces our results as far as the NUGM benchmark model is con- 
cerned. In this case, antiprotons will reach a discrimination sensitivity of 250 GeV, and 



24 



100 200 300 400 500 600 

P ' , ' , ' , ' , ' 1 10"' 

Pamela ly Funncl - 



Pamela 3y \~-v 



10 

n 10"" 7^ 



_q c 

10 



Antiprotons 

~ Positions ~ ~~^^-~ ^ 10 

lO' : ^ \ ^ \ ^ \ ^ \ ^ j 10" 

Visibility Line 

_o 

0^ ^^ r Direct Detection (op, SI) 

>> Neutrino Tlsc. (\i f. Sun) 

r^ , „-2 Antideuterons (Satellite Around Earth) 

en 

> 

\0'' r 



-10 ' 



lO"*" ' • ' • ' • ' • ' • J 

100 200 300 400 500 600 

Neutralino Mass (GeV) 

Figure 14: Future exclusions limits on tlie Funnel model. Symbols and conventions are the same 
as in fig. |l^. 

positrons of 200 GeV, after 3 years of Pamela data taking. For comparison, we see that 
antideuteron searches will do as good as antiprotons, while direct dark matter experiments 
will probably be marginally able to exclude models up to 400 GeV, though the visibility 
ratio lies in this case so close to 1 that this does not guarantee that XENON would see any 
signal, even at very low masses. 

6. The halo model dependence 

All results presented so far are in the context of a rather conservative halo profile. The 
spherically symmetric Burkert model we have implemented, from the point of view of 
Galactic dynamics, gives a subdominant term in a Galaxy embedding a stellar disc close 
to maximal, i.e. a dark term which, in the inner portion of the Galaxy, is hardly providing 
any dynamical effect at all. From the point of view of structure formation, we have already 
mentioned that some mechanism has to be invoked to reconcile this final configuration 
with a large core with models with large cusp describing, according to the simulations, the 
Galaxy before the baryon infall. 

Of course, this is not the only consistent picture. In the regime in which the baryons 
settle in the inner portion of the Galaxy through a smooth and slow process, actually the 
back-reaction on the dark matter particle is expected to go in the opposite direction with 
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Figure 15: Future exclusions limits on the NUGM model. Symbols and conventions are the same 
as in fig. |l^. 

respect to what we have assumed so far. This is the limit of adiabatic contraction of the 
system, with no net transfer of angular momentum between baryonic and non-baryonic 
terms; assuming spherical density profiles and unchanged local velocity distribution, the 
mass distributions in the initial and final configurations are related by [p^: 



Mi{ri)ri = [Mb{rf) + McDM{rf)] rj 



(6.1) 



Here, Mi{r), Mi,(r) and Mcdm{i") refer, respectively, to the mass profile of the halo before 
the baryon infall (i.e. the form one can infer from N-body simulation results), the baryon 
component as observed in the Galaxy today, and the cold dark matter component in its 
nowadays configuration, with the CDM cusp which is preserved or, actually, increased. 
It is possible to derive models for the Galaxy within this framework and fully consistent 
with dynamical measurements, see, e.g., []69| , |4^ . We consider here, as an example, the 
adiabatically contracted profile derived implementing the CDM profile found in ref. |70| 
interpolating on the results from one of the simulations with highest resolution so far. As 
the profile is derived numerically we cannot give its explicit form here; we just mention 
that it is obtained for a halo profile with virial mass and concentration parameter of, 
respectively, Myir = 1.8 x 10^'^Mq and Cmr = 12, and that the local density of the final 
spherically symmetric dark halo is equal to PAro3(''o) = 0.38 GeV cm~^. The approximation 
of adiabatic contraction is assumed to be valid down to the radius of 1 pc, below which 
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Figure 16: A comparison of the current dark matter exclusion limits, on the AMSB model, for 
two different Halo Models, the Burkert profile (solid lines) and the Adiabatically Contracted 
profile (dashed lines). Right: direct (cr^_p, spin independent) and indirect (muons from neutralino 
annihilations in the center of the Sun) detection Visibility Ratios. Left: the x^^^ for Antiprotons 
and Positrons. 



it is assumed that process which as led to the formation of the central black hole in the 
Galaxy W^ has erased any eventual enhancement (i.e. we do not have a spike in the form 
predicted, e.g., in p2|). As for the previous halo model, this configuration as well has been 
found after implementing all available dynamical constraints, and by the corresponding 
velocity distribution is derived self- consistently (for details, see [pq] ). 



Fig. 16 shows how the results for the current exclusion limits on the AMSB bench- 



mark model are affected resorting to the two different halo models we have introduced. 
Direct dark matter detection and neutrino telescopes rates are largely unaffected in the 
adiabatically contracted profile, and the conclusions we drew for the Burkert profile still 
apply. Antimatter searches are instead largely boosted by the new profile, which features 
a larger dark matter matter densities towards the Galactic center: antiproton flux fits rule 
out, with current data models with neutralino masses as large as 470 GeV, and positrons 
up to 370 GeV. Turning to future perspectives, fig. |^ illustrates that, again, the situation 
is not dramatically different for direct detection and neutrino telescopes, as compared with 



what shown in fig. 13; on the other hand, while antideuterons would probe masses 100 
GeV larger than in the Burkert profile, strikingly enough we find that both antiproton and 
positron searches will probe, in even only 1 year of data taking, the whole parameter space 
at the PAMELA experiment! 

Since the adiabatically contracted profile is cuspy towards the galactic center, we expect 
that another indirect detection channel, namely the observation of gamma rays in the 
direction of the center of the Galaxy, could also probe a large fraction of the AMSB model 
parameter space. The EGRET experiment, on the Compton Gamma- Ray Observatory has 
resolved a 7-ray source towards the GC [Q, tentatively extended (~ 1.5°, of the order of 
the EGRET angular resolution) rather than point-like, and with a spectrum apparently 



27 



10 
10* 
10' 
10^ 

lO' 
.2 10° 

Bi 10 

1 10-^ 
g 10" 
10* 
10" 
10' 
10' 
10" 




Direct Detection, Burkert Profile 
Neutrino Tlsc, Burkert Profile 
Antitieuterons, Burkert Profile 
Direct Detection, Adiab. contr. Profile 
Neutrino Tlsc, Adiab. contr. Profile 
Antideuteruns, Adiab. cuntr. Profile 



200 300 400 

Neutralino Mass (GeV) 



10 
lO" 

10" 
10" 
10" 















V 


Antiprolons, Burkert Profile 






N - 


Positrons 


Burkert Profile 






s 


— Antiprotons, Adiabatically contracted Profile ^ 






- - Positrons 


Adiabatically contracted Profile 


N 






-.^ 


AMSB ; 




Pamela ly 


^ 


"^ 


--.. ~ --, 






^"~^^\^^^ 






Pamela 3y 




^~^^ 


~^^ 










^^^^ 



200 300 400 500 

Neutralino Mass (GeV) 



600 



Figure 17: Future dark matter discrimination perspectives on the AMSB model for the Burkert 
(solid lines) and for the Adiabatically Contracted (dashed hues) Halo Model profiles. The left 
panel shows the Visibility Ratios at XENON-1 ton direct detection and at ICECUBE, as in fig. |l^. 
We also indicate how discrimination perspectives from antideuteron searches on a satellite around 
Earth vary with the Halo Model. In the right panel we display instead the discrimination parameter 
Ifp for antiprotons and positrons, as well as the putative 1 and 3 years sensitivity of the PAMELA 
experiment (long dashed lines). 

incompatible with the spectrum expected for the diffuse 7-ray flux due to the interaction of 
primary cosmic rays with the interstellar medium, while possibly consistent with a WIMP- 



induced component |74|. Expected fluxes within our frameworks are compared with the 
intensity and spectrum of the EGRET 7-ray source; the reader should also keep in mind 
that alternative explanations for this source have been proposed |75|, as well as that it 
has been argued that the position of the source should not be identified with the Galactic 



center |76]. A clearer statement on both these points will be feasible with the upcoming 



measurements by the next gamma-ray mission in space, the GLAST satellite |77], and, 
eventually, it will be possible to derive even more stringent bounds. 



As shown in fig. 18, it turns out that the adiabatically contracted halo profile is greatly 



constrained by current data from gamma rays, which excludes, at the 2-a level, masses up 
to 600 GeV within the AMSB benchmark model, as well as puts significant constraints for 
the other frameworks. It should be noticed however that the signal in gamma-rays is totally 
dominated by contributions close to the Galactic center rather than evenly distributed along 
the line of sight, hence the signal is extremely sensitive to what has been assumed on the 
(essentially unconstrained) distribution of dark matter in the Galactic center region; in 
that respect, note that for the signal for the Burkert profile is suppressed by over four 
orders of magnitude. Antimatter fluxes are much less dependent on this specific feature, 
and halo model configurations implying an enhancement in the antimatter fluxes without 
overproducing gamma-rays are certainly viable. 

Finally, there are further effects inducing large enhancements for antimatter fluxes and 
affecting less critically the other signals. E.g., one can consider halo models with are flat- 
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Figure 18: The bound coming from the continuous gamma rays background as measured by 
EGRET in the bin with the largest sampled energies (4 GeV < i?^ < 10 GeV) , for the two halo 
models considered here and for the three benchmark scenarios. As a function of the Neutralino 
mass, we plot the ratio of the expected signal in the considered energy bin over the 2-(t upper bound 
from the actual measured datum. Models above the horizontal green line are not consistent with 
EGRET and the particular considered halo model, and are therefore ruled out. The upper lines 
represent the results for the adiabatically contracted halo profile, while the three lower lines those 
for the Burkert profile. 



tened towards the Galactic plane, rather than being spherical: an increase in the local halo 
density enters linearly in direct detection and neutrino telescope rates, vi^hile quadratically 
for cosmic ray fluxes (the effect of flattening has been sketched, e.g., in |^9|). Moreover, 
the dark matter halo may not be perfectly smooth but have a clumped component, again 
with possibly large enhancements in antimatter yields, see, e.g., |[7^. 



7. Conclusions 

We considered in this paper three benchmark scenarios, respectively featuring a bino, wino 
and a higgsino-like lightest neutralino. All models have large neutralino pair-annihilation 
cross sections, and the resulting relic abundance is compatible with the WMAP range 
provided some cosmological relic density enhancement mechanism, or non thermal neu- 
tralino production, is operating. We described in detail the spectral features of antiprotons, 
positrons and antideuterons for each of these models. 
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Making use of a self-consistent halo model, we first directly compared the current 
performance of different search strategies, namely direct spin- independent, neutrino flux 
from the Sun and indirect detection through positron and antiproton fluxes; we concluded 
that, in all cases, the only constraints presently coming from dark matter searches on the 
viable parameter space of the models under consideration stem from antimatter searches. 
In many instances, and although available experimental data on antimatter have a poor 
statistics, the efficiency of antimatter searches currently does better than direct detection 
by various orders of magnitude. 

Regarding perspectives for dark matter searches at forthcoming experimental facilities, 
we introduced a new parameter I^ which allows, given a SUSY model, to reliably assess 
its visibility at future experiments, without any bias from the kind of data analysis which 
will be carried out. In particular, we provided the expected critical values of /^ for the 
discrimination capability of the PAMELA experiment after one and three years of data 
taking. 

The comparison of future experimental dark matter search strategies in the three 
benchmark models considered once again shows that antimatter searches may be highly 
competitive with respect to both direct detection and neutrino telescopes. In some cases, 
such as for antiprotons in a wino dark matter scenario or for antideuterons within a 
mSUGRA funnel scenario, antimatter searches may be the only viable dark matter de- 
tection technique. 

Finally, we showed how our conclusions would be affected taking into account a more 
cuspy profile: while direct detection and neutrino telescope rates are hardly affected, anti- 
matter fluxes are largely enhanced. The typical increase in the value of the /^ parameter 
is between two to three orders of magnitude. Although the particular halo profile we con- 
sidered features large gamma rays yields, different models for the Dark Matter distribution 
at the very center of the Galaxy would suppress this signal, while the expected antimatter 
fluxes would realistically lie between the two extreme cases we discussed here. 

In the context of SUSY models with large annihilation cross sections, and in view of the 
imminent launch of space-based dedicated experiments, antimatter searches are therefore 
to be considered as a highly promising path towards the detection of dark matter. 
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